Study of intrinsic spin and orbital Hall effects in Pt 
based on a (6s, 6p, 5d) tight-binding model 
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We study the origin of the intrinsic spin Hall conductivity (SHC) and the d-orbital Hall conduc- 
tivity (OHC) in Pt based on a multiorbital tight-binding model with spin-orbit interaction. We 
find that the SHC exceeds fOOO he^-Q' 1 cm when the resistivity p is smaller than ~ 10 pit cm, 
whereas it decreases to 300 he -1 ■ Sl _1 cm _1 when p ~ 100 pQ, cm. In addition, the OHC is still 
larger than the SHC. The origin of the huge SHE and OHE in Pt is the large "effective magnetic 
flux" that is induced by the interorbital transition between d xy - and d ;c 2_ H 2-orbitals with the aid of 
the strong spin-orbit interaction. 
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Recently, the spin Hall effect (SHE) has attracted 
much attention due to its fundamental interest and its 
potential application in spintronics. The SHE has a close 
relation to the anomalous Hall effect (AHE) in ferromag- 
nets: In 1954, Karplus and Luttinger (KL) [1] studied 
the Hall effect in multiband systems and found that an 
electric field induces a spin-dependent transverse current 
in the presence of spin-orbit (SO) interaction. This effect 
causes the AHE (transverse charge current) in ferromag- 
netic metals and the SHE (transverse spin current) in 
paramagnetic metals. These phenomena are fundamen- 
tal issues in recent condensed matter physics [2-17]. In 
these years, great progress on the SHE in semiconductors 
has been made. Murakami et al. [2] and Sinova et al. [3] 
have studied the intrinsic (impurity-independent) SHE 
in semiconductors by developing the theory of KL. Now, 
the SHE in two-dimensional electron gas (2DEG) with a 
Rashba-type SO interaction is well understood [4-7]. Al- 
though the SHE in semiconductors was recognized by the 
optical detection of spin accumulation [8, 9], it is unfor- 
tunately too small for quantitative analysis. Therefore, 
materials that show a large SHE are highly desirable. 

Recent experiments have revealed that the SHE also 
exists in metals such as Al [10] and Cl-doped ZnSe [11]. 
In particular, the huge spin Hall conductivity (SHC) in 
Pt at room temperature [240 he^ 1 ■il _1 cm _1 ] [13], which 
is 10 4 times larger than the SHC reported in semiconduc- 
tors, has attracted great attention. Simple 2DEG mod- 
els cannot explain this experimental fact. Recently, the 
present authors have studied the SHE in Sr2Ru04 that 
is described by the i2 9 -orbital tight-binding model [14] 
and found that the anomalous velocity due to interorbital 
hopping gives rise to huge SHC in transition metals. This 
mechanism also causes the large AHE [15-17]. To reveal 
the origin of the huge SHE in Pt, we have to investigate 
the anomalous velocity due to the multiorbital effect by 
considering all the d-orbitals (t2 g +e g orbitals). 

In this letter, we study the intrinsic SHE and the d- 



orbital Hall effect (OHE) in Pt by analyzing a realistic 
multiorbital tight-binding model. In the low-resistivity 
regime where p < 10 //ficm, both the SHC and or- 
bital Hall conductivity (OHC) are constant of order 
1000 ~ 3000 he^ 1 ■ Sl^cm -1 , whereas they are strongly 
suppressed in the high-resistivity regime where p 3> 
10 pQcm. The derived coherent- incoherent crossover is a 
universal property of intrinsic Hall effects [15, 17]. Both 
the SHE and OHE originate from a kind of Peierls phase 
factor due to the "effective magnetic flux" [14] that is 
induced by a combination of the angular dependence of 
d-orbital wave functions and SO interaction. In Pt, the 
dominant contribution to the SHE is given by the d xy - 
orbital (in t2 g ) and the d x 2 _ y 2 -orbital (in e g ). Therefore, 
both the t2 g 

to explain the huge SHE in Pt 



and e g -orbitals should be taken into account 



Pt has a face-centered cubic (FCC) structure with 
a = 3.9 A, and the nearest inter-atomic distance is 2.8 
A. In the present study, we use the Naval Research Lab- 
oratory tight-binding (NRL-TB) model [18, 19] to de- 
scribe the bandstructure in Pt. The NRL-TB model 
employs the scheme of the two-center, non-orthogonal 
Slater-Koster (SK) Hamiltonian [20]. The SK param- 
eters are represented with distance- and environment- 
dependent parameters that are determined so that the 
obtained total energy and the band structures agree well 
with those obtained by the first-principles calculations. 
Here, we take into account 6s, 6p, and bd orbitals (in 
total, nine) and hopping integrals up to the sixth near- 
est neighbor sites. The electron number per Pt atom is 
ten. The NRL-TB model uses non-orthogonal bases, but 
we neglect the overlap integrals between different sites 
in this study for simplicity [21]. This simplification sig- 
nificantly changes the dispersion of the s-band far from 
the Fermi level, while that of the band structure near the 
Fermi energy is little affected. In the presence of the SO 
interaction for 5d electrons Hso = A X^('' s )i; the total 
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FIG. 1: (Color online) (a) Band structure of a (6s, 6p, 5d) 
tight-binding model for Pt (A = 0.04). T = (0,0,0), 
X= (7r,0,0), W = (tt, tt/2, 0), L = (7r/2,7r/2,7r/2), and 
K= (37r/4, 37r/4,0). (b) a* y (k) for 7 = 0.002 along T-X. 

(c) a^ y (k.) on the plane fc z = 0.197T, which is plotted only 
for trl^k) > 25. The crosses represent the Fermi surfaces. 

(d) fc^-dependence of for k y — k z — 0.197r; the minimum 
bandsplitting near the Fermi level A is ~ 0.035. 
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where Hg is a 9 x 9 matrix given by the NRL-TB model. 
The matrix elements of I are given in ref. [22]. The 
bandstructure obtained for the (6s, 6p, 5d) tight-binding 
model with A = 0.04 Ry is shown in Fig. 1 (a), which is 
in good agreement with the result of a relativistic first- 
principles calculation [23, 24] near the Fermi level. Ac- 
cording to optical spectroscopy, A = 0.03 Ry for a 5d 
electron in Pt, and A = 0.013 Ry for a Ad electron in 
Pd [22]. Hereafter, we set the unit of energy Ry; 1 Ry 
= 13.6 eV. Based on the NRL-TB model, spin wave ex- 
citations and the electron self-energy corrections in the 
ferromagnetic Fe are studied using the random-phase ap- 
proximation [25]. 

The 18 x 18 matrix form of the retarded Green function 
is given by G R (k, u>) = (ui+fi — H + iT)^ 1 , where \x is the 
chemical potential and T is the imaginary part of the k- 
independent self-energy (damping rate) due to scattering 



by local impurities (or inelastic scattering by phonons). 
The charge current in the present model is 
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Here, jjj = —e^ 3 -, where — e is the electron charge and 
fi — x,y. In this case, the atomic SO interaction is not 
involved in the charge current since it is k-independcnt. 
Then, the s z -spin current = {J^,s z }/2 is expressed 
as 
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Here, we discuss the current vertex correction (CVC) 
due to the local impurity potentials in the Born ap- 
proximation, which is given by AJ^/ oc ^ k G A J^G R . 
When a is one of the p-orbitals and /3 is one of the 
(s, d)-orbitals, (H (k)) a ^ is an odd function with re- 
spect to k «-> — k, and therefore the (a, /3)-component 
of (d/dk^jG — GJ^G is an even function. Note that 
\p v ) — » —\pu) (y = x,y,z) under the parity transforma- 
tion. Thus, (A J^)a,p is finite only when either a or (3 
is a p-orbital. In Pt, however, we have verified that the 
CVC affects the SHE only slightly [less than 5%] since 
the 6p-level is 20 eV higher than the Fermi level /j, and 
the p-electron density of states (DOS) at [i is very small. 
For this reason, we disregard the CVC hereafter. 

According to the linear response theory [26], the SHC 
is given by a z xy = a zI y + a^J , where 
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Here, / and II represent the "Fermi surface term" and 
the "Fermi sea term", respectively. In the same way, 
the OHC of the Fermi surface term O zI y and that of 
the Fermi sea term 0*y are given by eqs. (4) and 
(5), respectively, by replacing with the Z z -orbital cur- 
rent = {Jjj z }/2. Because of the cubic symmetry 
°f Pt , vftt, = °x V ■ e nvS and 0£„ = O z xy ■ e^ s , where 
S — x,y,z and e^s is the antisymmetrized tensor 
with e xyz = 1. 

When T a p = jSap (constant 7 approximation), ui- 
integration in eq. (5) can be performed analytically as 
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shown in ref. [15]: Then, a zI1 = a zIIa + a zIIb , where 
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where I, m represent the band indices. E^ is the Zth 
eigenenergy of H measured from the chemical potential 
M! UlaH^Upm = E l k 6i m , where a, /? are the orbital 
indices and U is a k-dependent unitary matrix. (Jf ) m/ 
in eqs. (6) and (7) is given by £ Q/3 t// Q ( J^U fjm . 

In the Born approximation, T in the Green function is 
given by n imp I 2 ± £k lm(G A (0) -G R (0)), where / is the 
local impurity potential and ni mp is the impurity density. 
When X/Ep <C 1, T is almost diagonal with respect to the 
orbital index; — "f a ^a,p [14]. In the case of Sr2Ru04, 
the SHC in the Born approximation is nearly three times 
greater than that in the constant 7 approximation [14]. 
In Pt, in contrast, we have verified that both approxi- 
mations give a similar SHC in the clean limit. For this 
reason, we use the constant 7 approximation hereafter. 

Here, we determine the part of the Fermi surface from 
which the SHC originates: In Fig. 1 (b) and (c), we show 

a£(k) ^ ^ <(k'), where <(k) 
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is the integrand in eq. (4). [Apparently, J2k a xyO i ) = 
&xy-\ &xy(k) is finite only on the Fermi surface, and 
it takes huge values at (0.737T, 0,0) (on T-X) and at 
(0.42-7T, 0.42-7T, 0.427r) (on L~r) since two bands are very 
close on the Fermi level in the present model. However, 
the contribution of these two points to the SHC is small 
after taking the k-summation. The dominant contribu- 
tion comes from a wide area around (0.19-7T, 0.197T, 0.577r) 
as shown in Fig. 1 (d). Here, the bandsplitting A near 
the Fermi level is 0.035. 

Now, we perform the numerical calculation of the SHC, 
using 128 3 ~ 512 3 k-meshes. Figure 2 (a) shows the 
A-dependence of the total SHC a z y = <r zI y + <j z y for 
7 = 0.02, which is smaller than A = 0.035. a zI1 rep- 
resents the Fermi sea term in eq. (5). Apparently, 
a %y ~ a xy ^> (Txy ls realized, a* increases with A mono- 
tonically, and it reaches 1000 her 1 -Sl^cm -1 at A = 0.03. 
To clarify the origin of the SHE, we study the SHC when 
the SO interaction is anisotropic: As shown in Fig. 2 
(a), a z xy for Hso = ^J2$xS x + lySy)i is much smaller 
than that in the isotropic case where Hso — A J2i(l' s )i]- 
On the other hand, a z y for Hso = A^i(^ s z)i almost 
coincides with that in the isotropic case. Therefore, it 
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FIG. 2: (Color online) (a) A-dependence of the SHC. The re- 
lation a% y ~ a^y 2> Uxy is realized. The matrix element of 
the SO interaction between the d^-orbital and the d x 2. y 2- 
orbital gives the dominant contribution to the SHC. Note 
that 1 [\e\/2na] « 1000 he' 1 ■ n^cm' 1 for a = 4 A. (b) 
Crossover behaviors of the SHC and OHC at 7 ~ A ~ 0.035. 
p is approximately proportional to 7. p = 1 corresponds to 
1000 pflcm. 



is concluded that the z-component of the SO interaction 
gives the decisive contribution to the SHC. The matrix el- 
ement of l z is finite only for (yz\l z \zx) = ~{zx\l z \yz) = i 
and (xy\l z \x 2 -y 2 ) = — (x 2 -y 2 \l z \xy) — 2i. Among them, 
the d xy - and <i a .2_ !/ 2-orbitals, both of which are given by 
the linear combinations of l z = ±2, cause a dominant 
contribution to the SHC as shown in Fig. 2 (a). 

Here, we discuss the 7-dependence of the SHC and 
OHC: When 7 is sufficiently small, these intrinsic Hall 
conductivities are proportional to the lifetime of the in- 
tcrband particle-hole excitation: H/A [1-3, 14, 15, 17]. In 
fact, Fig. 2 (b) shows that both the SHC and OHC for 
A = 0.03 [22] are independent of 7 for 7 < A - 0.035. In 
the high-resistivity regime where 7 > A, both SHC and 
OHC decrease drastically with 7 since the interband ex- 
citation is suppressed when the quasiparticlc lifetime h/j 
is shorter than hj A. This coherent-incoherent crossover 
of the intrinsic Hall conductivities {cr X y = const, for 
7 <C A and a xy oc p~ 2 for 7 ^> A) has been analyzed 
theoretically in refs. [15, 17]. In Pt, the SHC decreases 
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much faster than p~ 2 in the high-resistivity regime and 
the SHC becomes negative for 7 > 1, which may be due 
to a complex multiband structure. If we put 7 ~ 0.07, 
p ~ 0.1 [|e|/27ra] ~ 100 pflcm. Then, the obtained SHC 
is ~ 300 he^ 1 ■f2 _1 cm _1 , which is close to the experimen- 
tal SHC of Pt [13]. In the experimental situation, 7 in 
Fig. 2 (b) corresponds to h/2r within the spin diffusion 
length in Pt (~10 nm) from the interface of the junction, 
which might be larger than the bulk value of 7. 

We comment on <r xy b , which is frequently called the 
"Berry curvature term" . When T al 3 = -fSap and 7 — ► 0, 
a xy = a xy b ~ <7 xl [15]. However, a xy b is totally different 
from a xy in the high-resistivity regime as shown in Fig. 2 
(b), since the cancellation between a xy and <J xy a becomes 
worse when 7 is large. In many systems including Pt, 
a xy « a xy is realized for a wide range of parameters [15]. 

We briefly discuss the SHC using the Born approx- 
imation, where 7 Q is proportional to the DOS for the 
a-orbital, p a (0). When j a is a-dependent, a xy b ^ a xy 
even in the clean limit [14]. In fact, the SHC in Sr2Ru04 
given by the Born approximation is much larger than 
that given by the constant 7 approximation since the a- 
dependence of p a (0) is large [14]. In Pt, however, both 
approximations give similar results. For this reason, we 
use the constant 7 approximation. 
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FIG. 3: (Color online) Effective magnetic flux for ^-electron 
in the two-dimensional Pt model. This is the origin of the 
huge SHC and AHC in Pt. 

Figure 3 shows the FCC crystal structure of Pt on 
the xy-plane. Based on this two-dimensional model, we 
explain an intuitive reason why the huge SHC appears 
in Pt, by considering only d xy -, d x i. y 2-, and s-orbitals. 
±t represents the hopping integrals between the nearest 
neighbor d^-orbital and s-orbital, and ±t' is used for the 
next nearest neighbor d x s. y 2 -orbital and s-orbital. Both 
the hopping integrals change their signs by rotation by 
7r/2. Here, we consider the motion of a "f-spin electron 
on the left side of Fig. 3 along a triangle of a half unit 
cell: An electron in the d^y-orbital can transfer to the 
d x 2_ y 2 -orbital and vise versa using the SO interaction for 
a "[-electron h\l z /2; (xy\l z \x 2 -y 2 ) = — {x 2 -y 2 \l z \xy) = 2i. 



By considering the sign of the interorbital hopping in- 
tegral (±t and ±i') and matrix elements of the SO in- 
teraction, we can verify that a clockwise (anticlockwise) 
motion along any triangle path with the SO interac- 
tion causes the factor +i (— i). This factor can be in- 
terpreted as the Aharonov-Bohm phase factor e 27 ™^ 1 ^° 
[(/>o = /ic/|e|], where <j> represents the "effective magnetic 
flux" [14] <f) = § Arfr = ±0 o /4. This effective magnetic 
flux gives rise to the SHC of order 0{\). 

We also discuss the origin of the OHE by considering 
the motion of an electron with \l z = +2) oc \x 2 -y 2 ) + 
i\xy). We can show that an electron with \l z = ±2) in 
Pt acquires the Aharonov-Bohm phase, which gives rise 
to the OHC of order O(A ) [14]. 

In summary, we have studied the origin of huge SHC 
and OHC in Pt using a (6s, 6p, 5d) tight-binding model, 
and found that the SHC reaches 1000 he^ 1 ■ f^cm -1 
in the low-resistivity regime where p < 10 pQcm. Other 
significant findings of the present study are that (i) the 
OHC is still larger than the SHC in Pt, which will cause 
large surface magnetization of Pt; (ii) the huge SHC and 
OHC originate from the effective magnetic flux created 
by the d xy - and d a .2_ J( 2-orbitals; and (hi) the coherent- 
incoherent crossover behaviors of the SHC and OHC are 
derived by taking both the /-term and //-term into ac- 
count correctly. When p ~ 100 pflcm, the obtained 
SHC becomes comparable with the experimental value 
240 he^ 1 ■ f2 _1 cm — 1 . Note that p in the present calcula- 
tion corresponds to the resistivity within the spin diffu- 
sion length (~10 nm) from the interface of the junction. 
We comment that the effect of the overlap inregral re- 
duces the magnitude of the SHC for Pt to some extent 
[21]. Finally, we discuss the role of the Coulomb inter- 
action: Although the SHC is independent of the renor- 
malization factor z = (1 — <9£(w)/<9w) _1 | w= o (= m/m*) 
[14, 15], it will depend on the w-dependence of 7(0;) as 
well as the CVC due to the Coulomb interaction. They 
are important future issues. 

During the preparation of this paper, we found a pa- 
per where the SHC was calculated based on a relativistic 
first-principles calculation [27]. Only the Berry curva- 
ture term given by eq. (7) was calculated, which is not 
justified in the high-resistivity regime (7 3> A) [15]. 
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